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SUMMARY 

While sliding mode observers (SMOs) using discontinuous relays are widely analyzed, most SMOs 
are implemented computationally using a continuous approximation of the discontinuous relays. This 
approximation results in the formation of a boundary layer in a neighborhood of the sliding manifold in 
the observer error space. Therefore, it becomes necessary to develop methods for attenuating the effect 
of the boundary layer and guaranteeing performance bounds on the resulting state estimation error. In 
this paper, a method is proposed for constructing boundary-layer SMOs (BL-SMOs) with prescribed state 
estimation error bounds. The BL-SMO formulation is then extended to simultaneously estimate exogenous 
inputs (disturbance signals in the state and output vector fields), along with the system state. Two numerical 
examples are presented to illustrate the effectiveness of the proposed approach. Copyright (c) 2016 John 
Wiley & Sons, Ltd. 
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1. INTRODUCTION 

Estimation of system states and exogenous inputs (disturbance inputs in the state and output 
vector fields) for nonlinear systems is a critical problem in many applications. These applications 
include: estimating actuator faults in mechanical systems, unmodeled disturbances in biomedical 
systems, and attacks in the measurement channels of cyberphysical systems [1-3]. The presence 
of exogenous inputs generally degrade closed-loop system performance. Therefore, it is imperative 
to design observers that simultaneously estimate state and exogenous inputs for implementation of 
high-performance closed-loop control systems. 
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The application of sliding modes [4-9] to state and unknown input observer design has been 
widely developed in the context of linear systems. In [10] and [11], an equivalent output error 
injection term is proposed to recover the state and measurement disturbance signals. Linear matrix 
inequalities for the construction of the observer gains and the reconstruction of the state disturbances 
are proposed in [12-15] for linear systems. An extension to Lipschitz nonlinear systems has been 
proposed in [16-23], and one-sided Lipschitz nonlinear systems in [24]. These observers are only 
applicable for a limited class of nonlinear systems and furthermore, most of the above papers deal 
with the estimation of states and unknown inputs acting in the state vector field, in the absence of 
measurement disturbances. 

Descriptor systems provide an attractive approach for simultaneous estimation of the states 
and exogenous disturbances [25, 26]. Sliding mode observer based on descriptor systems is 
proposed in [27]. Some recent papers [28, 29] also discuss reconstruction of the unknown 
signals using second-order sliding modes. However, the classes of nonlinearities considered in the 
current literature are restricted to Lipschitz or quasi-Lipschitz nonlinearities, which may introduce 
conservativeness in the design. Additionally, the SMOs are formulated with a discontinuous 
injection term, but are implemented with a continuous (generally sigmoidal) injection term. This 
is because the system response with the discontinuous injection term is computationally taxing, and 
difficult to implement in practice. 

We ameliorate some of these open problems in the present paper. Our contributions include the 
following: (i) we propose a method to simultaneously reconstruct the system state and exogenous 
inputs (both disturbances in the state and output vector fields); (ii) we provide ultimate bounds on 
the observer estimation error based off a tractable continuous approximation of the discontinuous 
relay term (this is sometimes called a ‘boundary-layer’, see for example [30-33]); (iii) we extend 
current SMO formulations to a wider class of nonlinearities using incremental multiplier matrices, 
proposed in [34]; and, (iv) we demonstrate the utility of smooth window functions in recovering the 
exogenous inputs to within a prescribed accuracy. 

The rest of the paper is organized as follows. In Section 2, we provide our notation. In Section 3, 
we define the class of nonlinear systems considered and formally state the objective of this paper. 
Subsequently, an observer architecture is presented and sufficient conditions are provided which, if 
satisfied, specify performance bounds on the observation error of the plant state and unknown output 
disturbance signal. In Section 4, we leverage smooth window functions to reconstruct the unknown 
state disturbance signal to within a prescribed accuracy. In Section 5, we test our proposed observer 
scheme on two numerical examples, and offer conclusions in Section 6. 


2. NOTATION 


We denote by M the set of real numbers, and N denotes the set of natural numbers. Let p G N. 


For a function / : M i-a M, we denote the space of p-times differentiable functions. The function 


/ G if 
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and / G Coo if sup]^ |/| < oo. For every v G we denote \\v\\ = Vv^v, where is the transpose 
of V. The sup-norm or oo-norm is defined as ||^||oo — sup^^j^ ||i;(t)||. We denote by Amin(F') 
the smallest eigenvalue of a square matrix P. The symbol P {-<) indicates positive (negative) 
definiteness and Ay B implies A — B y 0 for A, B of appropriate dimensions. Similarly, bb) 
implies positive (negative) semi-definiteness. The operator norm is denoted ||P|| and is defined as 
the maximum singular value of P. For a symmetric matrix, we use the ^ notation to imply symmetric 
terms, that is, 

a b 
c 



a b 


1 


For Lebesgue integrable functions g,h,wo use the symbol * to denote the convolution operator, that 
is, 

/ oo /‘OO 

h{t — r)g{t) dr = h{t)g{t — r) dr. 

-oo J —oo 


3. PROBLEM STATEMENT AND PROPOSED SOLUTION 
We begin by describing the class of systems considered in the paper. 


3.1. Plant model and problem statement 
We consider a nonlinear plant modeled by 

x{t) = Ax{t) + Bff{t, u, y, q) -|- Bgg{t, u, y) + (la) 

y{t) = Cx{t) + Dwy{t). (lb) 

Here, x = x(t) G is the state vector, u = u{t) G is the control action vector, y = y(t) G 
'MAy is the vector of measured outputs. The nonlinear function g = g{t,u,y) :'K x x 'MAy \-y 
models nonlinearities in the system whose arguments are available at each time instant t. 

Let the function / = /(t, u,y,q) x 'MA'^ x 'MAy x 'MAy \-y denote the system nonlineari¬ 
ties whose argument q G R^« has to be estimated, where 

q = CqX, 


and Cq G R’^'^ x R’^^^ . 

The signal Wx — Wx{t) G R"^^ is the unknown state disturbance, for example: unmodeled 
dynamics, actuator faults or attack vectors. The signal Wy = Wy(t) G R"^^ models unknown 
measurement/sensor disturbances, for example: cyber-attacks on the measurement channel or 
sensor faults. We refer to the vectors Wx and Wy as the exogenous input. The matrices A, Bg, Bf, 
G, G and D are of appropriate dimensions. 

To proceed, we make the following assumptions. 

Assumption 1 

The right-hand-side of (la) is locally Lipschitz. 
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Assumption 2 

The matrices G and D have full column rank, that is, rank(G) = rrix and rank(I)) = rriy. 
Assumption 3 

The state disturbance input Wx{t) is Lebesgue integrable. 

Finally, we make an assumption on the classes of nonlinearities considered in this paper. To this 
end, we need the following definition. 


Definition 1 (Incremental Multiplier Matrix) 

A matrix M G is an incremental multiplier matrix if it satisfies an incremental 

quadratic constraint ((^QC) 

r n T r n 


6q 

Sf 



> 0 , 


( 2 ) 


where 

Sq^q^-q^e 


(3a) 


and 

Sf = /(t, u, y, qi) - f{t, u, y, ^ 2 ) G (3b) 


for all (t, u, y, gi, ^ 2 ) G M x x W^y x x 


To illustrate the concept of the incremental quadratic constraint, we provide the following 
examples. 


Example 1 

Consider the nonlinearity f{t,y,q) = cos q. Since 


I cosgi - cosg 2 | < \qi - ^ 2 !, 


we have 


(gi -q 2 f - (cos gi - cos^ 2 )^ > 0 , 


that is. 


r -n 

T 



p -l 

1 

to 


1 0 


1 

to 

cos gi — cos q 2 


0 -1 


cos gi — cos q 2 


Hence, an incremental multiplier matrix for / is 


M = 



0 

-1 


( 4 ) 


for any C > 0- 
Example 2 

Consider the nonlinearity f{t,y,q) = g|g|, which is not globally Lipschitz. The nonlinearity / 
satisfies the inequality 

(91I91I - 921^2 I) (91 - 92) > 0 , 


Copyright ©2016 John Wiley & Sons, Ltd. 
Prepared using rncauth.cis 


Int. J. Robust. Nonlinear Control (2016) 
DOI: 10.1002/mc 














STATE AND EXOGENOUS INPUT RECONSTRUCTION USING BL-SMO 


5 


for any qi,q 2 C M. This can be rewritten as 


qi - q2 

T 

0 

1 

qi - q2 

_qi\qi \ - q2\q2\_ 


1 

0 

_qi\qi \ - q2\q2\_ 


Hence, an incremental multiplier matrix for f{q) is 


M = 



1 

0 


(5) 


for any C > 0- 
Remark 1 

Clearly, if a nonlinearity has a non-zero incremental multiplier matrix, it is not unique. Any positive 
scalar multiplier of an incremental multiplier matrix is also an incremental multiplier matrix. 

The class of nonlinearities satisfying ^QC contains a wide class of nonlinearities, including 
globally- and one-sided Lipschitz nonlinearities, incrementally sector-bounded nonlinearities, 
incrementally positively real nonlinearities, and nonlinearities with derivatives lying in cones or 
poly topes. For more details regarding the construction of the incremental multiplier matrix for 
different categories of nonlinearities, we refer the reader to the Appendix. For the more interested 
reader, we refer to [34,35]. 


Unknown 

exogenous 

inputs Wy-, Wy 



Figure 1. Overview of proposed state and exogenous input estimation scheme. The unknown exogenous 
inputs Wx , Wy acting on the nonlinear plant are shown in red. The estimated state x and estimated exogenous 

inputs Wy,Wx are shown in blue. 


3.2. Objective and overview of the proposed solution 

Our objective is to design a boundary-layer sliding mode observer (BL-SMO) for the nonlinear 
system (1) that can simultaneously reconstruct the state x{t), and the unknown exogenous inputs 
Wx{t),Wy{t). We first rewrite the nonlinear plant (1) as a generalized descriptor system whose state 
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is the augmented state 


n T 


w' 


Subsequently, we use the descriptor system as a platform to 


design a BL-SMO and guarantee ultimate bounds on the estimation error of the augmented state, 
that is, the estimation error of x and Wy. We conclude this section by demonstrating that the BL- 
SMO error dynamics converge to the boundary-layer sliding manifold in finite time, which will be 
used in the sequel to reconstruct the state disturbance input Wx. For convenience, an overview of the 
scheme is provided in Figure 1 . 


3.3. Generalized descriptor formulation 


Let 


- A 
X = 


Wy 


be an augmented state vector. Also let 

E = 

A = 


Iflx 0 

A 0 


C= C D e 




^Ux X inx-\-my) 


"bTlx X (tT-x -{-TTly ) 


^nyX{nx+my) 


Then we can represent the nonlinear plant (1) as a descriptor system 


Ex{t) = Ax{t) + Bff{t, u, y, CqEx) + Bgg{t, u, y) + Gw^it), 
y{t) = Cx{t). 


We illustrate this with an example. 


(6a) 

(6b) 


Example 3 

Suppose we have the nonlinear system 

xi u 

Xi + 3X2 — X2-\-Wx 

Clearly, we can write this in the form (1) with 


, y = xiEWy. 


A = 


1 0 
1 3 


Bf = 


0 

-1 


5 Bg — 


,G = 

Thus, the descriptor system can be written as 


C = 


r n T 
1 


,D = fCg = 


r n T 
0 


f = xl, g = u. 


1 

0 

0 


X 


1 

0 

0 


X 


1 


0 

o 

0 




X 






= 






+ 


u + 


xl + 


II 

1 0 1 



0 

1 

0 


_Wy_ 


1 

3 

0 


Wy_ 


0 


-1 


1 



Wy_ 


□ 
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Remark 2 

Note that the descriptor system (6) is constructed for developing the observer. It is not 
computationally implemented. For simulations, we use the original nonlinear plant (1). 

To proceed, we require the following technical result. 


Lemma 1 

Suppose the number of measured outputs is greater than or equal to the number of sensor 
disturbances; that is, riy > rriy. Then there exist two matrices Ti G and T2 G 

j^(na,+my)xny SUCh that 

TiE-T2C = In^+my- (7) 


Proof 
LetT = 


Ti T 2 


and 


E 


' I o' 

-C 


-C -D 


Computing T reduces to solving the linear equation TV = /. By Assumption 2, we know that D 
has full column rank, which implies V has full column rank. Hence, a left inverse of V exists. We 
denote as a left inverse of V, that is, V^V = /. Clearly, T = 1/^ is a solution to TV = /. 

Therefore, Ti can be computed by taking the first Ux columns of V^ and T 2 is the matrix 
constructed using the last Uy columns of V^. This concludes the proof. □ 


Remark 3 

A particular choice of such a left inverse is the Moore-Penrose pseudoinverse, that is, V^ = 

{V^V)-^V^. □ 


3.4. Proposed BL-SMO 
Let 

ey = y - Cx. 

We propose the following boundary layer sliding mode observer architecture to estimate the plant 
states X and the exogeneous inputs Wx and Wy\ 

z = QzT{Li- QT2)y + T^Byg + T^Bff + TiGw^ 
x = z-T 2 y 

f = f{t, u, y, CqEx + L2ey) 

pFey/\\Fey\\ if \\Fey\\>T] 

pFcy/ri if IIFe^ll < ?7, 

where y is an available (measured) output, is a continuous injection term for the sliding mode 
observer parametrized by the smoothing coefficent 7 ^ > 0 and 

Q = TiA- LiC. 

The signal w^{t) will be used in the sequel to recover the state disturbance input Wx{t). 



(8a) 

(8b) 

(8c) 

(8d) 
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The observer is parameterized by four gain terms: (i) the linear gain Li G (h) the 

innovation term L 2 G which improves the estimate of the known nonlinearity / by adding 

a degree of freedom in the design methodology, (hi) the matrix F G and, (iv) the scalar 

p > 0. 

Remark 4 

Assumption 1 implies that the observer ODEs also have unique classical solutions as G and 
hence, the functions x, x are absolutely continuous. □ 


3.5. Derivation of error dynamics 

We investigate the error dynamics of the proposed observer. To this end, we first require the 
following result which is easily proven by verification. 

Lemma 2 

Let Q = TiA — LiC and R = Li — QT 2 , where Ti, T 2 are constructed as described in Lemma 1. 
Then TiA - QT 2 C - RC - Q = 0. □ 

We define the observer error to be 

e = X — X. 

Using (7) and (8), the observer error dynamics are given by 

e = X — X 
= X — Z T2CX 
= TiEx — z 

= Qe + {TiA - QT 2 C - RC - Q)x + TiBf{f - /) + TiG{w^ - 

Using Lemma 2 yields 

^ = (TiA - LiC)e + TiBfif - f) + TiG(u;, - wf). (9) 

Our objective is to design the observer gains Li, L 2 and F to ensure that the error dynamical 
system (9) is ultimately bounded and the effect of the unknown input Wx is attenuated. 


3.6. Ultimate boundedness of observer error dynamics 

In order to investigate the stability properties of the observer error (9), we need the following 
technical lemma. 

Lemma 3 

Suppose M = is an incremental multiplier matrix (see Definition 1) for the nonlinearity / and 
let 
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where / is defined in (8c). Then the condition 


CqE-L2C 0 
0 I 


M 


CqE-L2C 0 
0 I 


e >0 


holds for any x, i G . 


Proof 

Recall that y = Cx. From (6a) and (8c), we have 


f - f = f{t, u, y, CqEx) - fit, u, y, CqEx + L 2 {y - Cx)). 

Let qi = CqEx, q 2 = CqEx + L 2 iy — Cx) = CqEx + I/2(7(x — x), Sq = qi — q 2 , and 6f = 
f{t, u, y,qi) — f{t, u, y,q 2 )’ Hence, we obtain 6q = {CqE — L 2 C)e. Now, we can write 


5q 


'CqE-L2C o' 

e 


'CqE-L 2 C 0 ’ 

N. 


- 1 

0 

_1 

N. 


1 

0 

1_ 




( 10 ) 


Recalling that the matrix M is an incremental multiplier matrix of /, and substituting (10) into the 
incremental quadratic constraint (2), we obtain the desired matrix inequality. □ 


Herein, we present sufficient conditions in the form of matrix inequalities for the observer design. 
Theorem 1 

Let \\wxi')\\oo ^ Px and a > 0. Suppose there exist matrices Li, L 2 , E, P = y 0, an 
incremental multiplier matrix M = for the nonlinearity /, and scalars p, /i > 0, which satisfy 


S H- ^ ^ 0, 

G^T^P = EC, 
P I 


I y 


P 0 

P > Px^ 


(lla) 

(llb) 

(llc) 

(lid) 


where 


P - C^Yi + PTiA - YiC + 2aP PTiBf 
BjT^P 0 


and $ 


CqE-L2C 0 


0 / 

then the error trajectories of the BL-SMO (8) with gains Li = P~^Yi, L 2 , E, and p satisfies 


limsup ||e(t)|| < 
t^oo V ^ 


PVPx 


( 12 ) 


where p > 0 is the smoothing coefficient of the continuous injection term w'^. 
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Proof 

With Yi = PLi, we get 

_ _ \{TiA - Licyp + P{TiA - LiC) + 2aP * 
BjT^P 0 

Now, we consider a quadratic function of the form 

vm) = eitfpm- 


Herein, for readability, we omit the argument of e{t). Then, the time derivative of V(e) evaluated on 
the trajectories of the error dynamical system (9) is given by 


V(e) = 2e'^P{TiA - LiC)e + 2e^PTiBf{f - f) + 2e'^PTiG{w^ - w^). 


Let^ = 



if-fV 


. Then from (11a), we get 


0 + 

= 2e'^P{TiA - LiC)e + 2e'^PTiBf{f - f) + 2ae^Pe + 
= V{e) + 2aV{e) + M^i - 2e^PTiG{w^ - w^). 

From Lemma 3, we know that > 0. Hence, 

V(e) + 2aV(e) - 2e'^PTiG'(«;^ - w^) < 0. 


For error states satisfying ||FCe|| > rj, we have 

V(e) < -2Q;V(e) + 2e^PTiG{w^ - w^) 

< -2aV(e) + 2||u;^||||e"^PTiG|| - 2e^PTiGw^. (13) 


Hence, recalling the definition of from (8d) and condition (11b), we get 

F{y - Ci) 


V(e) < -2aV(e) + 2pA\e~^PTiG\\ - 2pe^PTiG 
= -2aV(e) + 2pA\e^PTiG\\ - 2pe^PT^G 
= -2aV{e) + 2pA\e'^PTiG\\ - 2p 


\\Fiy-Cx)\\ 

FCe 


\\FCe\\ 
\\G~^T^Pe\\^ 


WG^T^PeW 

= -2aV{e) + 2pA\e'^PTiG\\ - 2p\\G^T^Pe\\ 
= -2aV{e) + 2||e^PriG||(p, - p). 


By choosing p > can ensure 


V(e) < -2aV{e), 


(14) 
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which implies global exponential stability of the observer error e to the set ||FCe|| < 77 with decay 
rate a; see for example, [31,36], for global exponential stability to a set. 

Now consider error states that satisfy ||FCe|| < 77 . Then, from (13), we obtain 

FCp 

V < -2aV{e) + 2\\G'^T^Pe\\\\w:,\\ - 2e'^PTiG—— 

= -2aV{e) + 2||G'’^r7^’e||||u;:,|| - 2e^PTiG ^^ 

V 

\\nTrpT Pp||2 

= -2aV(e) + 2\\G^T^Pe\\\\w,\\ - 2^-i-^ 

<-2aV{e) + 2\\G^T^Pe\\\\w,\\. 


Using (11b) yields 


V < - 2 aV(e) + 2 ||FC'e||||M;^|| < - 2 aV(e) + 2 ?? ||m;^||. 


Recalling that ||wj;|| < ||iea;(-)||oo < Px, we obtain 

V(e) < -2aV(e) + 2ripx. 


(15) 


Summarizing, we write 


V < 


-2aV 

—2cxV + 2ripx 


if ll-Feyll >7] 
if \\Fey\\ < p. 


The above implies that for any e G the inequality (15) holds. 

Note that taking Schur complements of (11c) yields / ^ pP. Hence e^e < pV(e). We use this 
inequality and the Bellman-Gronwall inequality on (15). This yields 


||e||2 < /re-2“(*-‘o)V(e(to)) + (1 - . 

2a \ J 


Taking the limit superior concludes the proof. 


□ 


Remark 5 

With L 2 and a fixed, the conditions in Theorem 1 devolve into a convex programming problem 
in P, Fi, P, M, /i and p. Additionally, solving the convex problem min^>o p subject to the 
constraints ( 11 ) with fixed P 2 and a results in tighter ultimate bounds on ||e||. 

Methods for converting the matrix inequality (11a) into LMIs without pre-fixing L 2 are provided 
in [35]. 


Remark 6 

Suppose the conditions of Theorem 1 are satisfied. Let x = 
the following holds for the plant state estimation error: 



X and w 




X. Then 


limsup \\x{t) - x{t)\\ < ^/pr]px/a, (16a) 

t^oo 
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and the measurement noise estimation error satisfies 

limsup \\wy{t) - Wy{t)\\ < 117]px/a. (16b) 

t—yoo 

For a fixed p, a and L 2 , we can minimize p over the space of feasible solutions. This 
attenuates the effect of p, thereby producing more accurate estimates of the state and measurement 
disturbance/sensor attack vectors. □ 

Remark 7 

Note that lim sup^^^ ||e|| 0 as 7 ^ -> 0, which implies that under ideal sliding (77 = 0) the matched 

disturbance can be completely rejected, and exact estimates of the plant state x and output 
disturbance Wy can be obtained. □ 

Summarizing, we have discussed a method to obtain estimates of the state x and measurement 
disturbance Wy to a specified degree of accuracy. However, certain applications such as fault 
detection [37] and attack detection [3, 38, 39], require the estimation of the unknown state 
disturbance input Wx> The following subsection provides a crucial ingredient for the simultaneous 
recovery of Wx along with Wy and x. 

3. 7. Finite time convergence to the boundary-layer sliding manifold 

We will demonstrate that the trajectories of the error system (9) are driven to the boundary layer 
sliding manifold in finite time. We begin with the following assumption on the plant states. 

Assumption 4 

The state vector x{t) and sensor disturbance Wy{t) of the nonlinear plant (1) are bounded, and 
known. 

The boundedness of plant states is reasonable for any practical system. We believe the restriction 
placed on the output disturbance is also not conservative, as measurement channels will transmit 
bounded signals, and attack vectors will be of finite magnitude. 

Herein, for brevity, we consider / to be a nonlinear function with the single argument q G . 
We use the following definition from [40, p. 406] to proceed with the development of technical 
results in this section. 

Definition 2 (Minimal Modulus of Continuity) 

The minimal modulus of continuity for any nonlinearity Lp{q) is given by 

7 y(r) = sup{||(^(gi) - (p{q 2 )\\ ■ 91,92 e ||9i - 92 II < r}, 


for all r > 0 . 

Remark 8 

An important property of the modulus of continuity is that it is a non-decreasing function, that is, if 
0 < ri < r 2 then 7 ^ (ri) < 70 (r 2 ). This follows from the definition of the supremum. □ 

We also pose a restriction on the class of nonlinearities considered. 

Assumption 5 

The nonlinearity f{q) considered in the plant (1) is uniformly continuous on 
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Remark 9 

Most nonlinearities encountered in practical applications adhere to Assumption 5. For example, if a 
function / is continuously differentiable with bounded derivative, Holder continuous with exponent 
f3 G (0,1], or globally Lipschitz continuous (which is a very common assumption in the literature), 
then / is also uniformly continuous. Hence, our assumption is not restrictive. □ 

We now present the following technical result. 

Lemma 4 

If f{q) is uniformly continuous on then 7 /(r) ^ 0 as r ^ 0. 


Proof 

Let 5 > 0. By uniform continuity, there exists > 0 such that \\qi — ^ 2 || ^ 5 forces ||/(gi) - 
/(^ 2 )|| < This implies sup{||/(gi) — f{q 2 )\\ • ||^i — q 2 \\ < ^ which, in turn, implies that 

^{6) < £. Since 7 (') is non-decreasing, we have 7 (r) < 5 for all r G [0,S], which concludes the 
proof. □ 


We are now ready to state and prove the following theorem which provides conditions for the 
observer error trajectories to converge to the boundary-layer sliding manifold 

5^ = {e G : ||FCe|| < r]} 


in finite time. 


Theorem 2 
Let 


Sf = fiCgEx) - f{C^M - L 2 {y - (71)), 


S = FC, a = Se, and 


Ai = A^in(G^T7PTiG). 


(17) 

(18) 


Suppose Assumptions 1-5 hold. If there exists matrices Yi, I/ 2 , F, P and scalars //, a which satisfy 
the conditions (11). If p is chosen to satisfy 


Ai/ 9 > sup IIS' {{TiA - LiC)e + TiGw^ + TiBjSf) ||, (19) 

t>to 

then the BL-SMO (8) with gains Pi = P~^Yi, P 2 , P and p generates error trajectories e{t) that 
converge to Sr^ in finite time. 


Proof 

If IlSell < 7/, we are done. Hence, for the remainder of this proof, we consider error trajectories 
satisfying UPeH > 7 ^. It is enough to show that 
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for some C > 0 in order to prove finite-time convergence to 5^, as argued in [41]. To this end, 


= a~^Se 

= ((TiA - LiC)e + TiBfSf + TjGK “ ^2)) 

Qp 

= a'^S {(TiA - LiC)e + TiBfSf + TiGw,) - pe^S'^STiG-— 

\\oe\\ 

Qp 

= a'^S {(TiA - LiC)e + TiBfSf + TiGw,) - pe^S'^(T iGYP{T iG)^^ 


from (11b). From [29], we know that for (11b) to have a solution, TiG must be full column rank. 
Hence 

(TiG)'^P(TiG) 0, (20) 

since P ^ 0. 

Recalling that Ai is the minimal eigenvalue of the symmetric positive definite matrix 
(TiG)Tp(TiG), we get 


a < S [{TiA 


LiC)e + TiBfSf + TiGwx) — p\i 


l|5e-|| 


< ||(t|| {\\S{{T^A-L^C)e + T^Bf5f + T^GwY\\-pXi)- 


We claim that for every C > 0, we can select p large enough to ensure a < —C||cr||. To prove our 
claim, we first demonstrate that 


sup ||PC {{TiA - LiC)e{t) + TiBfSf + TiGwYt)) || < oo. (21) 

t>to 

Using the triangle inequality, we have 
sup ||5 {{TiA - iiC)e(f) + TiGwYt) + T^BfSf) || 

t>to 

< sup (||P(TiA - Ti(7)||||e(t)|| + ||PTiG|| ||u;,(f)|| + \\ST,Bf\\\\f{C,Ex) - f{C,Ex - L 2 {y - (7l))||) 

t>to 

< sup \\S{T,A-L,C)\\\\m\\ + \\ST,G\\p^ + \\ST,Bf\\jf{\\CgE-L2CMm\\), 

t>to 

since Wx {t) and e{t) are bounded by Assumptions 3 and 4 and by (12). From (14), we also know that 
||e(t)|| decays exponentially when llPell >r]. This implies that ||e(t)|| is bounded and decreasing 
with increasing t. By Remark 8, this implies that Yf{\\CqE — L 2 C'||||e(t)||) also decreases with 
increasing t, since e(to) is bounded. Hence sup^>^Q 7 /(||CgP — P 2 C'|| ||e(t)||) is bounded. As all the 
terms are bounded, thereby finite, the condition (21) holds and the gain p selected using (19) is 
well-defined. □ 
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4. RECOVERING THE STATE DISTURBANCE USING SMOOTH WINDOW EUNCTIONS 


In this section, we discuss a filtering method to reconstruct the state disturbance Wx{t). Specifically, 
we will show that any piecewise uniformly continuous state disturbance input can be reconstructed 
to prescribed accuracy by filtering the injection term ( 8 d) of the BL-SMO with a smooth window 
function. 

To begin, we present our notion of a smooth window function. 

Definition 3 (Smooth Window Eunction) 

A smooth window function h{t) : R i-a M satisfies the following conditions: 

(i) h(t) G that is, h is smooth; 

(ii) h{t) > 0 for t G [—1,1] and h{t) = 0 elsewhere; 

(iii) = 

In the following proposition, we demonstrate that a uniformly continuous function can be 
approximated arbitrarily closely by filtering with smooth window functions. 

Proposition 1 

Let h{t) be a smooth window function, and define 




( 22 ) 


Also, suppose ti,T2 G R and : (n, T2) ^ R is uniformly continuous. Then for every 5 > 0, there 
exists dL P > 0 such that 

( 23 ) 

for every t G [n + /3, T 2 — /3]. Here, denotes the convolution operator. 


Proof 

We begin by noting that conditions (i)-(iii) in Definition 3 imply that the function / 1/3 is non-negative 
on the compact support [—/3, /3], and. 


hpit) 


(24) 


Using the definition of convolution, we have 


i’{t) - 


f 


i>{t - T)hp{T) 


dr. 


Applying (24) to the above yields 


V’(i) - 


f 




- T))hp{T)dT. 


Since by definition, is uniformly continuous on [n, r 2 ], for every 5 > 0 there exists a /3 > 0 such 
that for any tifi 2 C [n, r 2 ] satisfying |ti — ^ 2 ! < /3, we obtain — 'ip{t 2 )\ < £. 
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Therefore, we get the estimate 

/ c» 

-'lp{t-T))his{T)\dT 

-OO 

/ OO 

|■^(^) -'tp{t-T)\\hp{T)\dT 

-OO 

h/dir) dr, 


' —OO 
POO 


< S 


since is non-negative. Applying (24) on the above, we get (23). This concludes the proof. □ 


We wish to extend the result in Proposition 1 to a more general class of functions, namely 
piecewise uniformly continous function, defined next. 

Definition 4 (Piecewise uniformly continuous) 

A signal Wx is piecewise (or sectionally) uniformly continuous if 

(i) the signal Wx exhibits finite (in magnitude) jump discontinuities at abscissae of discontinuity 
denoted 

r ={ti:ie I}, 

where I = {1,2,...}. Specifically, I is the set of integers i satisfying a < i < bfor some a < b, 
where a may be — oc and b may be oo, and ti < whenever 1 G I. 

(ii) there exists a scalar c > 0 such that - ti\> c for every i G I; 

(iii) the signal Wx is uniformly continuous on the closure of each open interval in M \ I and this 
uniformity is independent of the interval. More formally, for every e > 0 there exists ^ > 0 
such that if ri < r2 are in M \ I satisfying |ri — T2 | < and such that there is no i G / with 
Ti <ti< T2, then \\wx{ti) - Wx{t2)\\ < e. 

Assumption 6 

The state disturbance input Wx is piecewise uniformly continuous. 

Remark 10 

Assumption 6 ensures that unknown input signals do not exhibit Zeno behavior (infinite number 
of jumps in finite time intervals), which is a reasonable assumption for state disturbances such as 
actuator faults or unmodeled inputs in physiological systems. □ 

We are now ready to extend Proposition 1 to piecewise uniformly continuous Wx . 

Proposition 2 

Let (/? : M i-A M satisfy Assumption 6 and p G Coo- Let the function hj^it) be defined as in (22) and 
let 

Ip = \J[ti-p,U+p] ( 25 ) 

iei 
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denote the union of closed neighborhoods around each abscissa of discontinuity of ^p. Then for every 
5 > 0, there exists a 0 < /3 < c such that 





for t G M \ X /3 
for t e Ip. 


As before, denotes the convolution operator and || • ||oo is the Coo norm. 


(26) 


Proof 

We begin by noting that the existence of the constant c in Definition 4 implies that T is a set of 
Lebesgue measure zero. Hence, the convolution integrals over R are well-defined. 

First, fix 5 > 0 and recall the definition of hp{t) in (22). Since the function hp has compact 
support [—/3, /3], the convolution integral is evaluated over the window of length 2/3. Thus, we can 
directly apply Proposition 1 with = ip\^t._pp.^^_^pj to obtain (23) for any interval in Ip. By (iii) 
in Definition 4, we can select /3 independent ofiel. From this, we conclude 

lk(i) - {hi}*ip){t)\\ < e, 


fort eR\Ip. 

However, the same cannot be said for the points t e Ip because the function cp is not uniformly 
continuous across the point of discontinuity. Since P < c,wc know that the function jumps just once 
in the interval (t^ — /3, t^ -f /3). Then we write 




{ip{t') - ip{t' -T))hp{T) dr 


< r ||^(t')-^(t'-r)||||t.^(r)||dr 
J-/3 


^ 2||(/::?||oo 


since hp{r) dr = 1 by construction. This concludes the proof. 


□ 


Proposition 2 implies that for piecewise uniformly continuous Wx, the filtering approach using 
smooth windows leads to high accuracy reconstructions in all but neighborhoods of the points of 
jump discontinuity. 

The following theorem is the main result of this section. It is an extension of a low-pass filtering 
method proposed in [42], which was for linear systems with sliding manifolds of co-dimension one. 


Theorem 3 

Suppose Assumptions 1-6 hold, and there exists a feasible solution (P,Li,I/ 2 ,F,M,a,/i) 
satisfying the conditions (11) in Theorem 1. Let p be selected as in (19) and Ip be defined as 
in (25). Then for a given 5 > 0, there exist scalars , Pn^ > 0, a sufficiently large T > 0, a 

sufficiently small 7 ^ > 0 and a low-pass filter 


hp{t) = diag 



Priw ^ 



(27) 
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such that 


for all t > T. 




2px + 5 / 2 , 


for t e Ip 

for t G [T, 00 ) \ Ip 


Proof 

We begin by fixing 5 > 0 and choosing > 0 for 1 < /c < such that hp. (t) satisfies (26) for the 
jth component of Wx. We define 


ILp = diag 


1 

dhfi^ 

1 


■ 

/3i 

dt 

1 ■ ■ • 

dt 

1. 


(28) 


where || • ||i denotes the £1 norm. 

Let xi = \\npl X 2 = \\FC{TiA-LiC)l and X 3 = \\FCTiBfl with || • || denoting the 
operator norm. We note that for a given y^/^’s and 5 , we can choose si sufficiently small to ensure 
that 

ma.x{xiei,X2£i,X3lfi\\CqE - L 2 C\\ei)} e 

- 1 - - a’ 

Ai 6 

where Ai = Xmin{G^Tf PTiG) as defined in (18). Recall that Ai > 0, as G^Tf PTiG y 0. Note 
that by construction PC = G'^Tf P, which implies 


Ai = X^UFCTiG). 


(30) 


Let S = FC, and ts be the time at which the error trajectories enter the boundary layer sliding 
manifold. We know that < 00 as p satisfies (19) in Theorem 2. Furthermore, we know that e(') is 
an absolutely continuous function (see Remark 4). 

Therefore, we can apply integration by parts for t > ts and use the compact support and 
smoothness of hp to obtain 


/ OO nOC 

hp{t — r)Se{r) dr = / hp{t — r)Se{r) dr, 

-00 J—00 


which implies 


ntrp 

/ hp{t — r)Se{r) dr = / hp(t — r)Se{r) dr. 

Jt-p Jt-p 

Let 6f be defined as in (17). Replacing the error-derivative ^ using (9) gives 


t+/3 


/ hp{t — r)Se{r) dr = 

Jt-p 

nt+P ^ ^ 

/ hp{t — r)S{TiA — LiC)e{r) dr 

Jt-p 

nt+P 

+ / hp{t-r)STiBf8fdr 

Jt-p 

pt+p 

+ / hffit - t)STiG {Wxir) - w^{t)) dr. 

Jt-B 


(31) 
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We now rewrite the last term in (31) as 


{his * STsGw^m - (his * STsGw^m 


nt^P ^ rip-p 

/ hi 3 {t — T)Se{r) dr — / h^{t — r)STiBf6f dr 

Jt—j3 


nt-\-(3 


r 

/ h^{t — r)S{TiA — LiC)e{r) dr. 
Jt-^ 


(32) 


From Theorem 1 , we know that 


limsup ||e(t)|| < 

t^OO 


and from Theorem 2 we get ||5'e|| < for t > Thus, for a given f3 > 0 and > 0 chosen as 
in (29), there exists a sufficiently small 77 > 0 , and sufficiently large T > ts for which 

sup ||5'e(r)|| < and sup ||e(r)||<5i (33) 

Te[t-(3,t+i3] Te[t-(3,t+(3] 


for all t > T. 

The inequality (33) along with Lemma 4 implies 
sup \\Sf\\ 

= sup \\f{CgEx{T)) - fiGgE^r) + L 2 Ce(r))|| 

re[t-/3,t-\-/3] 

< sup II 7 / {CqEe{T) - L2Ce{T)) || 

re[t—/3,t+/3] 

<^f{\\CgE-L2C\\es). 

Therefore, using (24) and (28), we upper bound the right hand side terms in (32). That is, 

/ oo 

h^{t — r)Se{r) dr 

-00 


< WUpW sup ||S'e(T) 

Te[t—/3,t+^] 

< Xiei, 


(34a) 


/ OO 

h(i{t - pSTiBf (/(g(r)) - f{q{T))) dr 

-OO 


hp{t — r)S{TiA — LiC)e{r) dr 


< sup \\STsBf\\\\f{q{r))-f{q{r))\\ 

Te[t-/3,t+/3] 

<X3Vi\\CgE-L2C\\eP, (34b) 

<\\S{T,A-LsC)\\ sup ||e(r)|| 

< X 2 ei, (34c) 


for t > T. We know that STiG is symmetric positive definite, and hence, Ai > 0. Therefore, 

\\ST,G{wPt)-w2{m 


\\Wx{t) - w2{t)\\ < 


Ai 


(35) 
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Applying to the above (32) and (34) produces 

\\{hi3*w^){t) - {hi}*w2){t)\\ < hpit) * \\w^{t) - wl{t)\\ 

< hpjt) * \\STiG{w^{t) - w2{t))\\ 

Ai 

^ (xi + X 2 )si + X3lf{\\CqE - L2C\\ei) 

Ai 


By construction of si in (29), we get 


\\{hp*w^){t) - {hp*w2){t)\\ < e/2. 


(36) 


Recall the definition of from (25). We now use (26) in Proposition 2 and (36) to obtain 


\\Woo{t) - {hi3^W^){t)\\ 

= \\wj,{t) - {h/3 * + {h/3 * - {hp * 

< \\wx{t) - {hp * 'u;^)(t)|| + \\{hp * Wx){t) - {hp * w'^){t)\\ 


< 


‘^Px H“ ^/2, 


for t eXf 3 

for t G [T, (X)) \ X /3 


for t > T > ts- and r] sufficiently small. This concludes the proof. 


□ 


Theorem 3 implies the existence a bank of smooth window filters capable of reconstructing the 
vector-valued signal Wx up to arbitrary accuracy in all but neighborhoods of jump discontinuities. 


5. SIMULATION RESULTS 

In this section, the performance of the proposed observer-based state and exogenous input estimation 
formalism is tested on two numerical examples. The first example is a practical system, where 
the nonlinearity is globally Lipschitz continuous and there is one state disturbance and one output 
disturbance signal. The second example is a randomly generated system (to demonstrate the non- 
conservativeness of our approach) with multiple exogenous inputs and a non-Lipschitz nonlinearity. 


5.7. Example 1 

We use the single joint flexible robot described in [29] to test our observer design methodology. The 
nonlinear plant is modeled as in (1) with system matrices 


0 

1 

0 

0 

-3.75 

-0.0015 

3.75 

0 

0 

0 

0 

1 

3.75 

0 

-3.75 

-0.0013 
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Bf — — 


0 

-1.1104 

0 

1 


. G = B, 


■ 1 ■ 

0.5 

, C = 

'l 

0 

0 

o' 


' 0 ' 

0 

0 

1 

0 

, D = 

1 

0 


_0 

0 

0 

1 


-2 

_l-3. 









Here, the nonlinearity / = cos(x 2 ) is globally Lipschitz and its argument, X 2 , is not measured 



Figure 2. Simulation Results. (Top left) The unmeasured variable X 2 (t) is shown in blue, and the dashed 
red line is the estimated trajectory X 2 {t). We note that the estimate is satisfactorily close to the actual. (Top 
right) The error e(t) is plotted in blue with the dashed black lines showing the error bound computed to be 
0.082. (Bottom Left) The estimate of the state disturbance input Wx{t) shown after 40 s. Note that the low 

pass filtered estimate is highly accurate. 


directly. The function g = 2.3sin(xi) is known at all t >to because xi is a measured output. 
The control input is set to zero. Hence rix = 4, = 3, rrix = 1, rriy = 1 and rif = 1 and Cq = 

0 10 0. Thus f{q) = cos{q) with q = CqX = X 2 . From (5), we deduce that this nonlinearity 

has an incremental multiplier matrix 


M = C 


1 0 

0 -1 


where C > 0* 
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We select a = 0.5 and L 2 = -16.55 -90.07 80.54 . Using CVX [43], we obtain a feasible 

solution^ to the LMIs in (11), namely C = 0.09 and // = 18.5268, the matrix 


the observer gain 


P = 


25.55 

-6.01 

-10.51 

-8.53 

-2.32 

-6.01 

1.61 

2.76 

2.26 

0.52 

-10.51 

2.76 

15.74 

7.67 

-1.54 

-8.53 

2.26 

7.67 

5.19 

-1.43 

-2.32 

0.52 

-1.54 

-1.43 

3.81 


1.58 

-0.43 

- 0 . 21 " 



29.81 

-75.99 

-37.99 


Ti = 

5.33 

-14.25 

-7.12 

7 


-20.95 

62.13 

31.06 



-9.29 

27.93 

13.97 



^ = [ 1.32 

0.21 

0.2 . 



and the sliding surface matrix 


For simulation purposes, we consider a randomly generated initial condition 


(to) = 2.09 -2.17 -0.31 -8.58 


and the exogenous inputs are chosen to be = sawtooth(2t + 1) and Wy = square(4t). Hence, 
px = ^- The observer is initialized at 2 ; = 0 and the sliding mode gain is set at p = 100. Finally, the 
boundary layer sliding mode injection term is computed using p = 10“^. From Theorem 1, we 
get the error state bound 

limsup ||e(t)|| < ,~ 0.073. 

t^oo V Of 

A 9th-order Butterworth low-pass filter with window length /3 = 0.24 s is used to 
obtain the actuator fault signal estimate. The corresponding MATLAB implementation is 
butter(9, 0.12, 'low' ). The simulation results are shown in Figure 2. We compute the 
experimental mean squared error \\wx{t) — w^{t)\\‘^ « 1.59 X 10 ^ fromt G [20,80], which verifies 
that our reconstruction is highly accurate. 


5.2. Example 2 

To demonstrate that our assumptions are not restrictive, we will test our method on a randomly 
generated system of the form (1) with multiple unknown inputs and a non-Lipschitz nonlinearity. 


tWe find that minimizing the norm of Yi, hence Li, usually enables faster runtimes using MATLAB’s odelSs or 
ode23s. 
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Here, 


2.44 

5.32 

9.29 

8.63 ■ 

1.1 

-4.11 

1.82 

2.53 

-0.09 

0.9 

-2.91 

0.06 

-4.53 

-3.45 

-8.59 

-12.14 


Bn = 0,Bu = 0, 


G = 


" 0.04 

1.77 ■ 


■ 0 ■ 













'1 

0 

0 

o' 


1 

1.37 

0.3 

, Bf = 

-1 

, c = 

0 

0 

1 

0 

, D = 

0 

-6.14 

-0.56 

0 


_0 

0 

0 

1 


-1 


_-2.71 

0.05 _ 


_ 1 _ 









and the nonlinearity is / = X 2 |x 2 |. We set = 0 2 0 0 . Since the nonlinearity is globally 

Lipschitz, we know from (4) that it is characterized by an incremental multiplier matrix of the form 


M = C 


0 1 
1 0 


for some C > 0- We fix a = 0.5, L 2 = —0.04 —0.23 1.42 and use CVX to obtain ( = 12.59, 

/i = 4.71. 


P = 


8.45 

-1.41 

-1.12 

4.58 

5.99 

-1.41 

8.75 

4.35 

-1.2 

0.07 

-1.12 

4.35 

20.8 

-3.93 

-1.38 

4.58 

-1.2 

-3.93 

13.95 

-3.4 

5.99 

0.07 

-1.38 

-3.4 

11.26 


' 43.58 

-0.26 

43.58 ' 



14.94 

-2.12 

14.94 


Li = 

-8.19 

1.77 

-8.19 



-22.09 

-3.2 

-22.09 



-32.1 

1.18 

-32.1 


F = 

-10.13 

3.98 

-48.81 

-3.22 

-10.16 

4.01 



and 

We generate a random initial condition 

x{to) = [-32.94 -31.38 -26.19 -68.^ 


and the exogenous inputs are chosen to be 




3 cos(t) 

5 sawtooth(4t) 
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10 


< 5-5 

I 5 


^rror 

Irror Bound 


20 

t 


40 



Figure 3. Simulation Results. (Top left) The actual (blue) and estimated (red dashed) trajectories of the 
unmeasured state X 2 {t) are shown. (Top right) Zoom-in of error trajectory and the computed plant state 
error bound. We note that the bound (black dashed) is not conservative. (Bottom) We illustrate that the 
exogenous inputs are estimated with high accuracy. 


and Wy = 10 sm{3t), respectively. Hence, ^ 5.831 and ||rr’^(')||cx) < 10; therefore, these 
exogenous inputs are of significant magnitude. The observer is initialized at z = 0 and the sliding 
mode gain is chosen p = 200. The continuous injection term w'^ is computed with p = 10“^. 
Therefore, from Theorem 1 , we get the error state bound 


limsup||e(f)|| < J « 0.1048. 

t^OO V ^ 

Two different smoothing filters with /3i = 0.3 s, /32 = 0.1 s are used to obtain estimates of the 
unknown state input Wx{t). The MATLAB command to implement these smooth window filters 
is: smooth (in jectionTerm, ' loess' ) ;. The simulation results are shown in Figure 3. 
We note that although the unknown exogenous inputs are reconstructed with high accuracy, the 
sawtooth input Wx 2 exhibits overshoots and undershoots at the points of jump discontinuities, as 
predicted by Theorem 3. 


6. CONCLUSIONS 

We developed a methodology for constructing implementable boundary layer sliding mode 
observers for a wide class of nonlinear systems using incremental multiplier matrix as a unifying tool 
for the design procedure. We formulate linear matrix inequalities which, if satisfied, can be used to 
construct the observer with pre-specified ultimate bounds on the reconstruction error of plant states 
and unknown output disturbances. We also demonstrate the utility of smooth window functions in 
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recovering the unknown state disturbance signal and provide an upper bound on the exogenous input 
estimation error for state disturbance inputs exhibiting jump discontinuities, which has not been 
investigated previously. The proposed methodology has a variety of applications including fault 
detection and reconstruction for mechanical systems, high confidence estimation in cyberphysical 
systems and secure communication. 
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A. INCREMENTAL MULTIPLIER MATRICES EOR COMMON NONLINEARITIES 


We present systematic methods for the computation of incremental multiplier matrices for a variety of 
nonlinearities analyzed in this paper and encountered in practical systems. We refer the reader to [34, Section 
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6] for a detailed discussion of methods used to compute incremental multiplier matrices and corresponding 
derivations of these matrices. 

We begin by recalling the definition of Sq and Sf given in (3). 


A.l. Incrementally sector bounded nonlinearities 
An incrementally sector bounded nonlinearity satisfies the inequality 


{MiiSq + MuS/yX{M2iSq + M22<5/) > 0, 


(37) 


for some fixed matrices Mu, M 12 , M 21 , M 22 and for all X e A, where Af is a set of matrices. After 
representing the nonlinearity in the form (37), the incremental quadratic constraint (IQC) in (2) is satisfied 
by choosing 



where, 


Ma = MuXM 21 + mJiXMu, 
Mh = m7iXM22 + mJiX'^Mi2, 

Me = MJ2XM22 + M22X'^Mi2. 


A.2. Incrementally positively real nonlinearities 

For a class of incrementally positively real nonlinearities, that is, nonlinearities satisfying 


SfxSq > 0 , 


the corresponding incremental multiplier matrix is given by 



with /^ > 0. 


A. 3. Globally Lipschitz nonlinearities 

For a globally Lipschitz nonlinearity that satisfies \\Sf\\ < Lf\\6q\\ for some Lf > 0, we write 


{Lf5q + 5ffiLfSq-6f)>0 


and inequality (37) is satisfied by choosing 



with hv > 0. 
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AA. Quasi-Lipschitz nonlinearities 

Another class of nonlinearities considered in this paper is the so-called ‘one-sided’ or ‘quasi’ Lipschitz 
nonlinearities that satisfy 

Q6f < R6q^ 


for some G M, Q G and R = R^ G . An incremental multiplier matrix for this class of 

nonlinearities is given by 


M = n 


■2ZfR 


-Q 

0 


with K > 0. 


A. 5. Nonlinearities with derivatives residing in a poly tope 

Suppose we have a nonlinearity / that satisfies 


dl 

dq 


e 0 , 


where 6 is a polytope with vertex matrices 0i,... ,0r.ln other words, 

df 


dq 


= Hx), 


where 0{x) = Ylk=i Xk^k, and Xk satisfies Xk>0 for all fc € {1,... ,r} and J2k=i Xk = 1- Then a 
corresponding incremental multiplier matrix 


M = 


Mil 

m72 


Mi2 

M 22 


satisfies the matrix inequalities 


M22 ^ 0 

Mil + ^12^/c + ^12 + M 220 k h 0 


(38) 


for all /c = 1,..., r. An example of this class of nonlinearity is f{q) = [sin(gi) cos(g 2 )] ^ whose derivative 
cos(gi) 0 


-sin(gi)_ 


which lies in a polytope © with vertices 


0i = -O 2 


1 

0 


0 

0 


, and Os 


-O 4 


0 0 
0 1 


Another example that falls into this category is the Takagi-Sugeno fuzzy model, proposed in [44]. 


A.6. Nonlinearities with derivatives residing in a cone 

Suppose we have a nonlinearity / that satisfies 


dq 


G 
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where O is a cone with vertex matrices cji,..., cjr. In other words, 


dl 

dq 


^(x), 


where <^{x) = Yl'k=iXk^k^ and Xk satisfies Xfe > 0 for all fc € {1,... ,r}. Then a corresponding 
incremental multiplier matrix of the form (38) satisfies the matrix inequalities 


M220k = 0 
Mi20k + ojM i 2 y 0 


for all /c = 1 ,..., r. An example of this class of nonlinearity is f{q) = [qi 


1 

0 


^ 2 . 


which lies in a cone O with vertices 


^ 2 / 5 ], whose derivative is 


UJl 


1 

0 


0 

0 


and (jU 2 


1 0 
0 1 
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